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Abstract —In this paper, for the first time, we analytically 
prove that the uplink (UL) inter-cell Interference in frequency 
division multiple access (FDMA) small cell networks (SCNs) 
can he well approximated hy a lognormal distrihution under a 
certain condition. The lognormal approximation is vital because it 
allows tractable network performance analysis with closed-form 
expressions. The derived condition, under which the lognormal 
approximation applies, does not pose particular requirements on 
the shapes/sizes of user equipment (UE) distribution areas as 
in previous works. Instead, our results show that if a path loss 
related random variable (RV) associated with the UE distribution 
area, has a low ratio of the 3rd absolute moment to the variance, 
the lognormal approximation will hold. Analytical and simulation 
results show that the derived condition can be readily satisfied in 
future dense/ultra-dense SCNs, indicating that our conclusions 
are very useful for network performance analysis of the 5th 
generation (5G) systems with more general cell deployment 
beyond the widely used Poisson deployment. 

1. Introduction 

Small cell networks (SCNs) have been identified as one 
of the key enabling technologies in the 5th generation (5G) 
networks Q . In order to gain a deep theoretical understanding 
of the implications that SCNs bring about, new and more 
powerful network performance analysis techniques are being 
developed. In this context, new performance analysis tools can 
be broadly classified into two large groups, i.e., macro-scopic 
analysis and micro-scopic analysis [2-8]. 

The macro-scopic analysis assumes that both user equip¬ 
ments (UEs) and small cell base stations (BSs) are randomly 
placed in the network, often following homogeneous Poisson 
distributions [2,3]. The micro-scopic analysis is usually con¬ 
ducted assuming that UEs are randomly dropped and the BSs 
are deterministically deployed, i.e., the BS positions in the 
considered cellular network are known. Generally speaking, 
the macro-scopic analysis investigates network performance on 
a high level [2,3], while the micro-scopic analysis gives more 
detailed results for specific networks [4-8]. Note that both 
analyses are related to each other. The average performance 
of micro-scopic analyses conducted over a large number of 
realizations of BS deployments should be equal to that of the 
macro-scopic analysis, provided that the examined realizations 
of the deterministic BS deployments follow the BS distribution 
assumed in the corresponding macro-scopic analysis. 

In this paper, we focus on the micro-scopic analysis, and 
in particular, we consider an uplink (UL) frequency division 
multiple access (EDMA) SCN. Note that the analysis of a UL 


EDMA system is more involved than its downlink counterpart 
due to the power control mechanisms used at UEs. For the 
UL micro-scopic analysis, existing works either use 

• Approach 1, which provides closed-form but complicated 
analytical results for a network with a small number of 
interfering cells, each cell with a regularly-shaped UE 
distribution area, e.g., a disk or a hexagon Q; 

• Approach 2, which makes a empirical conjecture on the 
UL interference distribution and on that basis derives ana¬ 
lytical results for a network with multiple interfering cells 
placed on a regularly-shaped lattice, e.g., a hexagonal 
lattice 0, 

• Approach 3, which conducts system-level simulations to 
directly obtain empirical results for a complex network 
with practical deployment of multiple cells placed on 
irregular locations 0. 0. ®- 

Obviously, Approach 1 and Approach 3 lack generality and 
analytical rigorousness, respectively. Regarding Approach 2, it 
has been a number of years since an important conjecture was 
proposed and extensively used in performance analysis, which 
stated that the UL inter-cell interference with disk-shaped UE 
distribution areas could be well approximated by a lognor¬ 
mal distribution in code division multiple access (CDMA) 
SCNs and in FDMA SCNs |j^. This conjecture is vital 
since it allows tractable network performance analysis with 
closed-form expressions. However, there are two intriguing 
questions regarding this conjecture: (i) Will the approximation 
still hold for non-disk-shaped UE distribution areas? (ii) Will 
it depend on the sizes of the UE distribution areas? In this 
paper, we aim to answer those two questions, and thus making 
significant contributions to constructing a formal tool for the 
UL micro-scopic analysis of network performance. 

Compared with the previous works [5,6] of the micro-scopic 
network performance analysis based on empirical studies, the 
contributions of this paper are as follows: 

1) Our work, for the first time, analytically proves the 
conjecture in Q, i.e., the UL inter-cell interference in 
FDMA SCNs can be well approximated by a lognormal 
distribution under a certain condition. 

2) The derived condition, under which the lognormal 
approximation applies, does not rely on particular 
shapes/sizes of UE distribution areas and can be readily 
satisfied in future dense/ultra-dense SCNs, indicating 


that conclusions derived using this framework are very 
general and useful for network performance analysis. 

3) Based on our work, we propose a new approach to 
fill an important theoretical gap in the existing micro¬ 
scopic analysis, which either assumes very simple BS 
deployments or relies on empirical results. Specifically, 
we directly investigate a complex network with prac¬ 
tical deployment of multiple cells placed on irregular 
locations. In order to do that we provide a theoretical 
framework based on the lognormal approximation of 
the UL interference distribution, supported by rigorous 
theoretical analysis, as well as the conditions under 
which the approximation should be valid. 

II. Network Scenario and System Model 

In this paper, we consider UL transmissions and assume 
that in one frequency resource block (RB) at a time slot, only 
one UE is scheduled by each small cell BS to perform an UL 
transmission, which is a reasonable assumption in line with 
the 4th generation (4G) networks, i.e., the UL single-carrier 
FDMA (SC-FDMA) system and the UL orthogonal FDMA 
(OFDMA) system in the 3rd Generation Partnership Project 
(3GPP) Long Term Evolution (LTE) networks Q and in the 
Worldwide Interoperability for Microwave Access (WiMAX) 
networks fig, respectively. 

Regarding the network scenario, we consider a SCN with 
multiple small cells operating on the same carrier frequency 
as shown in Fig. [T] In Fig. [T] a total of B small cells 
exist in the network, including one small cell of interest 
denoted by Ci and B — \ interfering small cells denoted by 
Ch, b G {2,..., B}. We consider a particular frequency RB, 
and denote by Kb the active UE associated with small cell Cb 
using that frequency RB. Moreover, we denote by Rb the UE 
distribution area of small cell Cb, in which its associated UEs 
are distributed. 



Fig. 1. A schematic model of the considered SCN. 

The distance (in km) from the BS of Cb to the BS of Ci, 
b G B} and the distance from UE Kb to the BS of Cm, 

b,m G B} are denoted by Db and dbm, respectively. 

In this paper, we consider a deterministic deployment of BSs, 
and thus {Db} is assumed to be known. However, the position 
of UE Kb is assumed to be randomly and uniformly distributed 


in Rb- Hence, dbm is a random variable (RV), the distribution 
of which cannot be readily expressed in analytical form due 
to the arbitrary shape of Rb- 

Based on the definition of dbm, the path loss (in dB) from 
UE Kb to the BS of Cm is modeled as 

Bhm = A -f odo^-^^dbrn, (1) 

where A and a are the path loss at the reference distance 
dbm = 1 and the path loss exponent, respectively. Note that 
Lbm is a RV due to the randomness of dbm- In practice, A 
and a are constants obtainable from field tests ID- 

The shadow fading (in dB) from UE Kb to the BS of 
Cm is denoted by Sbm {b,m G B}), and it is usually 

modeled as a zero-mean Gaussian RV because the linear 
value of Sbm is commonly assumed to follow a lognormal 
distribution HI- Hence, in this paper, we model Sbm as an 
independently identical distributed (i.i.d.) zero-mean Gaussian 
RV with variance crfijad, expressed as Sbm ~ A/” (O, 

The multi-path fading channel vector from UE Kb to the 
BS of Cm is denoted by h.bm G where we assume 

that each UE and each BS are equipped with 1 and N omni¬ 
directional antennas, respectively. All channel coefficients are 
assumed to experience uncorrelated flat Rayleigh fading, and 
they are modeled as i.i.d. zero-mean circularly symmetric 
complex Gaussian (ZMCSCG) RVs with unit variance. 

The normalized maximal ratio combining (MRC) reception 
filter denoted by f^ = is considered at the BS of Cm, 

where is the Hermitian transpose of hmm- According 

to d). the effective channel gain associated with the link 
from UE Kb to the BS of Cm {b, m G [1,..., B} , b m), 
defined as \imii-bm\^ and denoted by Hbm, follows an i.i.d. 
exponential distribution exp(l) with the mean equal to 1. 

The transmit power of UE Kb is denoted by Pb- In practice, 
Pb is usually subject to a semi-static power control (PC) 
mechanism, e.g., the fractional pathloss compensation (FPC) 
scheme GD- Specifically, based on this FPC scheme, Pb in 
dBm is modeled aG] 

Pb = Pq P V {Lbb + Sbb) , (2) 

where Pq is the power basis in dBm on the considered 
frequency RB, rj G (0,1] is the FPC factor, Lbb has been 
defined in |^, and Sbb A/” (O, 

For clarity, the defined RVs in our system model are 
summarized in Table |I] 

*A more practical assumption would be the correlated shadow fading |l3) , 
which constructs Sbm and Sjm J'? ^ {1, • • • ? B} , h ^ j) as correlated 

RVs, and the correlation coefficient should decrease with the increase of the 
distance from UE Kb to UE Kj. Such assumption of the correlated shadow 
fading will greatly complicate the analysis since it is difficult to characterize 
the distribution of the inter-UE distance. For the sake of tractability, in this 
paper, we assume i.i.d. shadow fading for the UE-to-BS links. 

^In practice, Pb is also constrained by the maximum value of the UL power 
and affected by the per-UE signal-to-interference-plus-noise ratio (SINR) 
target. The power constraint is a minor issue for UEs in SCNs since they 
are generally not power-limited due to the close proximity of the UEs and the 
associated SCN BSs |12| . The per-UE SINR target, on the other hand, will 
greatly complicate the analysis since it is difficult to model the distribution 
of the target SINRs. For the sake of tractability, in this paper, we model Pb 
as 1^, which is widely adopted in the literature [3,6-8,11]. 





Table I 

Definition of RVs. 


RV 

Description 

Distribution 

dbm 

The distance from Ki, to Cm 

related to Rb 

Lbm 

The path loss from Kb to Cm 

related to Rb 

Sbm 

The shadow fading from Kb to Cm 

i.i.d. M (0, crfhad) 

^bm 

The channel gain from Kb to Cm 

i.i.d. exp (1) 

Pb 

The UL transmission power of Kb 

related to Rb 


III. Analysis of UL Interference 

Based on the definition of RVs listed in Table |I] the UL 
received interference power (in dBm) from UE to the BS 
of Cl can be written as 

h = Pb - Lti - Sbi + 10 logio Hbi 

= Pq + ivLbb — Lbl) + iv^bb — Sbl) + lOlog^g Hbl 

= Po + L + S + 10\og^oHbi, (3) 


where 0 is plugged into step (a) of 0- Besides, L and S 

are defined as L = {rjLbb — Lti) and S — {rjSbb — Sbi), 
respectively. In particular, since Stb and Sbi {b G {2,..., U}) 
are i.i.d. zero-mean Gaussian RVs, it is easy to show that S' 
is also a Gaussian RV with a distribution JV cr|), where 


(4) 


I MS = 0 

Vs = (1 + M^)<^lhad ' 

From the definition of /{, in Q, the aggregated interference 
power (in mW) from all interfering UEs to the BS of Ci is 


= ^10io 


h 


(5) 


6=2 


In the following, we analyze the distribution of I™"' 
in three steps. First, we investigate the distribution of 
(S + lOlogiQ Hbl) shown in 0. Second, we approximate the 
distribution of Ib as a Gaussian distribution. Third, we show 
that the distribution of can be further approximated as a 
lognormal distribution. 


A. The Distribution of {S + lOlogig Hbi) in 0 

According to the product of a lognormal RV and an 
exponential RV can be well approximated by another lognor¬ 
mal RV. Therefore, in our case, the sum {S + 10 logig i/bi) 
can be well approximated by a Gaussian RV G, because lOio'S’ 
is a lognormal RV and Hbi ~ exp (1). The mean and variance 
of G can be respectively computed as 0 


MG = Ms + Moffset 


( 6 ) 


— Og -I- yoffset 

In 0^ Moffset and cioffset are suggested to be —2.5 and 
5.57, respectively. Note that the approximation is very accurate 
when (t| > 36 Such requirement is readily satisfied in 
practical SCNs, e.g., it is recommended in Jl^ that — 
100 and hence Ug > 100 > 36 because of (|4]b 

B. The Distribution of Ib in 0 

In the following, we will prove that under a certain con¬ 
dition, Ib can be well approximated by a Gaussian RV. First, 


based on the approximation discussed in Subsection III-A 
reform 0 into 


we 


Note that the probability density function (PDF) and cumula¬ 
tive density function (CDF) of L are generally not tractable 
because L is a RV with respect to dbb and dbi, which jointly 
depend on the arbitrary shape of Rb- 

Next, we analyze the distribution of Ib by investigating 
the condition under which the sum of a Gaussian RV and 
an arbitrary RV, i.e., {L + G), can be well approximated by 
another Gaussian RV. To that end, we denote by /ig and 
the mean and variance of L, respectively, define a zero-mean 
RV as L = L — p,]^, and further define another zero-mean RV 
as G = G — PC- As a result, 0 can be re-written as 


Ib ^ L + G + (Pq + pl + Mg) • 


1) The distribution of ( L + G] in {8): 


( 8 ) 


of{l 

Obviously, if -f G^ 


a Gaussian RV, then Ib can a’ 


in ( 81 can be well approximated by 


so be well approximated by the 


same Gaussian RV with an offset (Pg + Ml + Mg)- In order 
to prove that 
a Gaussian R' 

Theorem [T] 


L + Gj indeed can be well approximated by 
we introduce the Berry-Esseen theorem in 


Theorem 1. [Theorem 1, 15] The Berry-Esseen theorem: Let 
Xi, X 2 , ■ ■ ■, Xm be M independent RVs with E{Aif} = 0, 
E = af > 0, and E < 00 , i G {1,2,..., M}. 

M 

Also, let X = y ' =. Denote by Fx (x) and T* (x) the 
CDF of X and the CDF of the standard normal distribution. 
Then, we have 

sup \Fx (x) - $ (x)| < Cofo, (9) 

icGK. 


where ipQ = 




and Gg = 0.56. 


In Theorem Gqi/jq is the upper bound of the maximum 
gap between Fx (x) and <I> (x). Therefore, Theorem shows 
how good the sum of M independent RVs can be approxi¬ 
mated by a Gaussian RV. As long as GgV'o is reasonably small, 
we can say that the approximation is good. Based on Theo¬ 
rem [T] we propose Femma S to address the question: Under 

in (8 I be well approximated by 


what condition can 
a Gaussian RV? 


(l + Gj i 


Lemma 2. 


Considering the zero-mean RV (^L -f cj given 
by {8), we define a normalized RV asY = Denote by 

V<yi+<yQ 

Fyyy) and $ (j/) the CDF ofY and the CDF of the standard 
normal distribution. Then, we have 

sup |Fy (j/) - (?/)| < T, (10) 


where r = 


Goe{|Z|''} 

{s/Pl+PgY 


and Gg = 0.56. 


Proof: In order to make use of Theorem [T] to prove 




























Lemma 1^ we construct M — 1 zero-mean Gaussian RVs by 
breaking the zero-mean Gaussian RV G into the sum of M — 1 
i.i.d. zero-mean Gaussian RVs {Gi} ,i G {1, 2,..., M — 1}, 
i.e., G = The variance of each G^ is cr^. = . 

Hence, the constructed M —I zero-mean Gaussian RVs {G^} 
and the RV L form the set of M i.i.d. RVs required to invoke 
Theorem [^ pejiote the normalized sum of the M RVs as 


Y = 


'f 


Go'00 Go X 


E 


Ef=7'E{|G.|^} 


(v^ 


E 


= GoX- 


L 


. ^^M—1 2 

+ Ei=l 


(M-l)(v^)^r(2) 


M-1 


where E 




(v^)''r(2) 



GnE 


Our proof is thus completed. 


An important note on the proof of Lemma is that the de¬ 
composition of G into the sum of M—1 i.i.d. zero-mean Gaus¬ 
sian RVs {Gi} allows the diminishing of E 

which quickly reduces GqV'o in Q to T, thus making it much 
easier for the Gaussian approximation to be valid. According 
to Theorem [T] if t in ( [T^ takes a reasonably small value, then 
(jL + G^ can be well approximated by a Gaussian RV. 

From Lemma l^^^we can see that r decreases with the 
decrease of E 


and with the increase of Gq 


^ and 


Therefore, if the 




approximation holds. Intuitively speaking, if 




.. Then, we can evaluate the upper bound 


= l “Gi 

of the maximum gap between Fy (y) and $ (y) by checking 
the metric GqV'o using Theorem [T] 

According to Theorem [T] the metric GqV'o for Y becomes 

3' 


fading. It is important to note that although the PDF and 
CDF of L and L are difficult to obtain, if not impossible, 

the values of (t|, ®||t| | can be easily computed using 

numerical integration over the UE distribution area Ri,. We 
will briefly discuss the calculation of r in the next subsection. 

To sum up, according to Lemma if t is reasonably small 
for the considered SCN, e.g., around 0.01 (an error about 1 
percentile), (^L -f G^ can be well approximated by a zero- 
mean Gaussian RV. Then, R in ([^ can be well approximated 
by a Gaussian RV Qb, the mean and the variance of which 
can be computed by 


— Pq + + MG 


(13) 


2 ) The calculation of t: 
Considering the definition of L, 


L respectively in (1^ 


and we can directly evaluate t using based on the 


-,( 11 ) 


has been 


results of and E( L ) computed by 


ML = 


/ 

JRt 


= / LdZ 


plugged into step (a) of ([n}. : Since the value of M is 
arbitrary in Theorem we let M approach infinity, and 
thus the first term of the numerator of i/'o in ([TT), i.e., 

(M-l)(C2)"r(2)/ 

IV M-i 

result Go-tpo converges to r given by 


= / {rjLbb — Lbi) dZ 
JRb 

dP 

{y- 1) A-hQlogio^ 


will diminish to zero and as a 


( 12 ) 


=/ 

JRb 

^1= [ {L- ycf 

JRb 

H 


dZ 


dV, 


(p- 1) A-halogio^ -pL 


dZ. 


dZ. 


E( L 


= / |t-ml| 

Rb 


dZ 






{y-l)A + alogio - fiL 

Rbl 


3 

dZ. 


(14) 


(15) 


(16) 


C. The Distribution of 7™'^ in (|^ 

According to GZ)’ the sum of multiple independent log¬ 
normal RVs can be well approximated by a lognormal RV. 
Considering the expression of 7™"' in if each Ib,b G 
can be approximated by the Gaussian RV Qb, 
we can conclude that 7™"^ can be well approximated by a 
lognormal RV, denoted by 7™"' = 10wQ. The RV Q is a 
Gaussian RV and its mean and the variance are denoted by 
fiQ and tTg, respectively. According to 117 , pg and erg are 


associated with UE distribution area 


obtained by solving the following equation set, 

B 

t'q Vi) = n ^Qb 


(si) = Gi 


Rij is low, then r tends to be a small value, and hence the 


6=2 

B 


t'q (s 2) — n 'i^Qb (*2) — G2 

6=2 


(17) 


is low. 


then the interfering UEs will be basically concentrated in a 
small area of Rb- As a result, the geometrical randomness 
of the interfering UEs is reduced, and hence the dB-scale 
shadow fading, usually modeled as Gaussian distributions, 
will dominate the dB-scale UL interference distribution. Also, 
a larger Gq allows a better approximation due to the more 
dominance of the Gaussian distribution of the dB-scale shadow 


where 'Fjc (s) is the approximated moment generating func¬ 
tion (MGE) evaluated at s for a lognormal RV defined as 
IOto^. Such approximated MGE is formulated as 


/ X Wm { {s/2G^am + h-X 

T'X [S) =2_^ -J= exp l-sexp|4G 

m—1 ^ \ \ 


C 


( 18 ) 



































where ( = Mg is the order of the Gauss-Hermite 

numerical integration, the weights Wm and abscissas Um for 
Mg up to 20 are tabulated in Table 25.10 in Usually, Mg 
is larger than 8 to achieve a good approximation GZl- 

In ([^, Si and S 2 are two design parameters for generating 
two equations that can determine the appropriate values of 
/ig and (Tg. Generally speaking, when si and S 2 take smaller 
(larger) values, the mismatch in the head (tail) portion between 
the actual CDF and the approximated one can be reduced GZl- 
The basic idea of GD is to let the two concerned RVs, 
i.e., and j™"', statistically match with each other in 

the sense of having two equal points on their MGFs. The 
solution of ([TtJ can be readily found by standard mathematical 
software programs such as MATLAB. 

Based on this lognormal approximation, the approximated 
PDF and CDF of can be respectively expressed as ||T^ 


f/mW (v) ! 




exp < - 


(Clnu - 


and 


^ 1 1 .,/C In f — Mo 

Fjmw (u) « - + -erf I -p= 

2^g 


(19) 


( 20 ) 


where erf (•) is the error function. 

D. Summary of the Proposed Analysis of UL Interference 
To sum up, we highlight the steps in our proposed analysis 
of UL interference in the following. First, we use ([T4]), Gg 
and (161 to calculate <j\ and e||^| | for ^ach Rh,b C 

{2,... ,i3}. Then, based on the calculated results, we check 
whether r given by Lemma is reasonably small. If it is the 
case, we can approximate /{, as a Gaussian RV Qfa according 
to the discussion in Subsection |III-B Finally, we approximate 
the RV I™"' in (|^ as a lognormal RV I™"' = lO^Q, with the 
distribution parameters ^g and Cg obtained from solving the 
equation set ( [l^ . The chain of approximation to obtain J™"' 
is summarized in Table [II] shown on the top of next page. 


IV. Simulation and Discussion 


In order to validate the results from the proposed micro¬ 
scopic analysis of UL interference, we conduct simulations 
considering two scenarios, with a single interfering cell and 
with multiple interfering cells, respectively. According to the 
3GPP standards GD’ the system parameters in our simulations 
are set to: A = 103.8, a = 20.9, Pg = —76 dBm, rj = 0.8, 
and cr| = 100. Besides, the minimum BS-to-UE distance is 
assumed to be 0.005 km G3- Regarding the equation set GZl' 
to determine /rg and ctq, we choose the parameters as si = 1, 
S 2 = 0.1 and Mg = 12 as recommended in fTT). 


A. The Scenario with a Single Interfering Cell 

In this scenario, we consider B — 2, as shown in Fig. 1^ The 
x-markers indicate BS locations, where the BS location of Ci 
has been explicitly pointed out. The dash-dot line indicate a 
reference disk to illustrate the reference size of small cell C 2 - 
The radius of such reference circle is denoted by r, and the 


distance between the BS of Ci and the BS of C 2 , i.e., D 2 , is 
assumed to be 1.5r. In our simulations, the values of r (in km) 
are set to 0.01, 0.02 and 0.04, respectively. In this scenario, the 
interfering UE K 2 is randomly and uniformly distributed in an 
irregularly shaped UE distribution area R 2 , as illustrated by 
the area outlined using solid line in Eig. The shape of i ?2 is 
the intersection of a square, a circle and an ellipse, which has 
a complicated generation function. Examples of the possible 
positions of K 2 within R 2 are shown as dots. 



Fig. 2. Illustration of R2 (irregular shape, r = 0.04). 


Table III 

Results of the proposed analysis (B = 2). 


The value of r 

T 

MQ 2 

^Q2 

0 

d 

II 

0.0082 

-97.10 

205.25 

r = 0.02 

0.0122 

-99.69 

207.73 

r = 0.04 

0.0184 

-101.54 

211.44 


Despite of the complicated shape of R 2 , our proposed 
lognormal approximation of the UL interference still works. 
Specifically, from ( G5| ), GH and Lemma we can 

calculate p.^, (t\, E< Z > and r, respectively. The results 


of r are tabulated for various values of r in Table m Erom 
Table |Inj we can observe that, when r = 0.01 or 0.02, the 
values of r are around 0.01. Consequently, it indicates a good 
approximation of I 2 as a Gaussian RV. Note that r = 0.01 
and r = 0.02 correspond to the typical network configurations 
for future dense/ultra-dense SCNs G)> which shows that the 
derived condition in Lemma can be readily used to study 
future dense/ultra-dense SCNs. 

To confirm the accuracy of the proposed approximation 
summarized in Table [Hj we plot the approximated analytical 
results and the simulation results of I 2 for the considered R 2 
in Eig. 1^ Also, the numerical results of and crg^ are 
shown in Table [111] As can be seen from Eig. the proposed 
Gaussian approximation of I 2 is indeed very tight for the 
considered i ?2 with such an irregular shape. Note that the 
approximation appears to be very good even for r = 0.04 
(despite a relatively large value of t at around 0.02), which 
implies that the derived condition in Lemma might be too 
strict and the proposed approximation might be extended to 
more general conditions. Hence, it is our future work to find 
a way to relax the sufficient condition proposed in Lemma 
to allow for a wider application of the proposed analysis. 


B. The Scenario with Multiple Interfering Cells 

The lognormal approximation of the UL interference for a 
network with multiple interfering cells placed on a hexagonal 


























Table II 

The chain of approximation to obtain 


Original RV 

Approximated RV 

Approximated distribution 

Distribution parameters 

S' -1- 10 login Hbi in i 

i 

G 

Gaussian 

See 0 

Ib ~ Po + -/j + G in 1 

7i 

Qb 

Gaussian 

See jiSi 

jmW ^ ill ,5, 

jmW 

Lognormal 

The solution of 1 17 1 



Fig. 3. The simulation and approximation of I 2 (irregular shape of R 2 ). 

lattice has been validated in 0- As a significant leap from 
the hexagonal network, in this subsection, we apply the pro¬ 
posed framework on a more complex network with practical 
deployment of multiple cells and provide the approximation of 
the UL interference distribution. 

Here, we consider a 3GPP-compliant scenario p^ , as 
shown in Fig. 0 In Fig- 0 B is set to 84 and all small cell 
BSs are represented by x-markers. Particularly, the BS of Ci 
has been explicitly pointed out. The reference UE distribution 
area for each small cell is a disk with a radius of r PH. In 
our simulations, the values of r (in km) are set to 0.01, 0.02 
and 0.04, respectively. The reference disk-shaped areas can be 
easily seen in Fig.0from any isolated small cell. However, due 
to the irregular positions of the cells, the actual UE distribution 
areas of the considered cells are of irregular shapes due to cell 
overlapping. The irregularly shaped UE distribution areas are 
outlined in Eig.0by solid lines. Interfering UEs are randomly 
and uniformly distributed in those areas. An important note 
is that the considered network scenario is different from that 
adopted in [2,3], where UE distribution areas are defined as 
Voronoi cells generated by the Poisson distributed BSs and 
those Voronoi cells cover the whole network area. In practice, 
small cells are mainly used for capacity boosting in specific 
populated areas, rather than provision of an umbrella coverage 
for all UEs. Therefore, the 3GPP standards recommend the 
hotspot SCN scenario depicted in Fig. 0 for UE distribution 
in the performance evaluation of practical SCNs. Nevertheless, 
the proposed micro-scopic analysis of UL interference can 
still be applied on a particular Voronoi tessellation, since the 
derived condition in Lemma 0 does not rely on particular 
shapes/sizes of UE distribution areas. 

To the best knowledge of the authors, the micro-scopic 
analysis of UL interference for such complex network as 
shown in Fig. 0 has never been attempted before in the 
literature. Now, we investigate the considered network with 
the proposed approach of performance analysis. First, we use 



Fig. 4. Illustration of a practical deployment of multiple cells (r = 0.04). 


14||, (|15|) and (|16|l to calculate p.^, (j\ and E 


I - |3 

\l\ 


for each 


Rb,b G {2, .. . 7 ^ 4 } displayed in Fig. 0 Then, based on the 
calculated results, we check the value of r for each Rb using 
Lemma 0 The maximum values of the 83 i?;,-specific r’s for 
various r values are presented in Table It is shown in 
the table that when r = 0.01 and r = 0.02 the maximum 
T is below or around 0.01, and thus each can be approx¬ 
imated by a Gaussian RV Qb according to the discussion in 
Subsection III-B Again, note that r = 0.01 and r = 0.02 
correspond to the typical network configurations for future 
dense/ultra-dense SCNs (TJ, which indicates the usefulness of 
the proposed approach of the micro-scopic analysis in future 
5G SCNs. Finally, we approximate the RV in (0 as a 
lognormal RV = lOraQ, with the distribution parameters 
PQ and tTg obtained from solving the equation set (117 1. The 
numerical results of pQ and Oq are provided in Table |lv 


are provided in Table IV 


Table IV 

Results of the proposed analysis {B = 84). 


The value of r 

Maximum r 

i^Q 


o 

d 

II 

0.0066 

-75.09 

17.77 

r = 0.02 

0.0125 

-77.21 

18.30 

r = 0.04 

0.0176 

-78.76 

18.55 


To validate the accuracy of our results on the UL inter¬ 
ference, we plot the approximated analytical results and the 
simulation results of in dBm for the considered network 
in Fig. 0 As can be seen from Fig. 0 the resulting lognormal 
approximation of is reasonably good for practical use. 
Note that the approximation shown in Fig. 0 is not as tight 
as that exhibited in Fig. 0 The noticeable approximation 
errors in Fig. 0 are mostly caused by the inaccuracy of 
approximating the sum of multiple lognormal RVs as a single 
lognormal RV in GD- Note that the parameters of si = 1 and 









































S 2 = 0.1 recommended by GZ) allow for an overall match 
between the actual CDF and the approximated lognormal CDF, 
though the approximation in the head portion and the tail 
portion of the CDF are compromised to some extent. Hence, 
a straightforward way to improve the approximation of UL 
interference shown in Fig.|^is to approximate the CDF of J™"' 
as a piece-wise lognormal CDFs, because smaller (larger) si 
and S 2 are helpful to increase the quality of approximation in 
the head (tail) portion of the considered CDF pT) . Regarding 
the individual approximation of If, as a Gaussian RV Qb, 
we find that it is tight for every considered Rb- Due to 
space limitation, we omit the investigation on the Gaussian 
approximation for each R, which is quite similar to the 
discussion in Subsection IIV-AI 



Fig. 5. The simulation and approximation of in dBm. 

A final note on the proposed network performance analysis 
is that our approach is very powerful to obtain the analytical 
results in an efficient manner. The computational complexity 
is mainly attributable to the numerical integration to obtain 

the values of and E||i| | for the 83 interfering 

cells. In contrast, the simulation involves a tremendously 
high complexity. Specifically, in our simulations, around one 
billion of realizations of R have been conducted for the 83 
interfering cells in order to go through the randomness of 
all the RVs listed in Table [I] This shows that the proposed 
micro-scopic analysis of network performance is elegant and 
computationally efficient, which makes it ideal to study future 
5G systems with general and dense small cell deployments. 

V. Conclusion 

The conjecture of approximating the UL inter-cell inter¬ 
ference in FDMA SCNs as a lognormal distribution is vital 
because it allows tractable network performance analysis with 
closed-form expressions. Compared with the previous works 
based on empirical studies, our work, for the first time, analyti¬ 
cally proved that the conjecture is conditionally correct and the 
derived condition does not rely on particular shapes/sizes of 
UE distribution areas. Based on our work, we proposed a new 
approach to directly and analytically investigate a complex 
network with practical deployment of multiple cells based on 
the approximation of the UL interference distribution. The 
proposed approach has the following merits. 


1) It can deal with any shape of UE distribution area, which 
tolerates more practical assumptions, e.g., irregular hot¬ 
spots, overlapped small cells, etc. 

2) It measures the quality of approximation using a closed- 
form expression and the simulation results validate the 
tightness of the approximation. 

3) It can cope with a large number of small cells by a low 
computational complexity of analysis. 

As future work, we will find a way to relax the derived 
condition of valid approximation to allow for a wider applica¬ 
tion of the proposed analysis. Besides, we will consider non- 
uniform UE distribution and other types of multi-path fading, 
e.g., Rician fading or Nakagami fading, in our analysis. Also, 
we will further investigate alternative distributions of the sum 
of multiple lognormal RVs to make our analysis more accurate. 
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